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This  study  addresses  the  effects  of  dynamic  strain  turnover  and  antibiotic  prophylaxis  on  rates  of  group  A  Streptococcus 
(GAS)  antibiotic  resistance  and  disease.  The  authors  analyzed  the  strain  distributions,  disease  rates  and  patterns  of 
antibiotic  resistance  of  802  GAS  isolates  collected  from  2002  through  2007.  These  samples  were  collected  from  patients 
with  GAS  infection  symptoms  at  ten  military  facilities.  Macrolide  resistance  peaked  at  25%  during  2004,  due  to  the 
geographically  widespread  dominance  of  a  single  resistant  strain  (M75).  The  resistant  strain  was  not  retained  regardless 
of  local  patterns  of  macrolide  use,  and  resistance  rates  decreased  upon  replacement  of  M75  with  macrolide-susceptible 
strains.  Disease  rates  were  similarly  correlated  with  dominance  of  specific  M  types.  Statistical  analysis  revealed  temporal 
correlations  between  strain  distributions  at  multiple  locations.  Only  the  most  common  strains  yielded  enough  data  at 
multiple  sites  for  statistically  significant  comparison  of  temporal  fluctuations  in  dominance,  but  these  (including  M44, 
M3,  Ml  8,  Ml  1 8  and  M6)  all  yielded  highly  significant  temporal  correlations  of  90%  or  greater  on  yearly  scales.  As  expected 
given  the  complexity  and  variability  of  strain  distributions  on  shorter  time  scales,  analysis  on  a  monthly  scale  yielded  lower 
degrees  of  positive  correlation  (31-62%),  but  in  this  case  all  significant  correlations  were  still  positive.  Shifts  in  antibiotic 
resistance  profiles  and  disease  rates  at  specific  sites  appear  to  be  associated  with  strain  replacements  happening  on 
larger  scales,  independent  of  antibiotic  use  at  individual  sites. 


Introduction 

Popularioris  in  crowded  conditions  for  extended  periods  of 
time,  exemplified  by  military  recruits  and  deployed  troops,  are 
extremely  susceptible  to  contagious  acute  respiratory  disease 
(ARD).*'^  The  agents  most  commonly  associated  with  ARD  in 
these  situations  include  group  A  Streptococcus  (GAS),  adenovi¬ 
rus  and  influenza.  Vaccines  provide  excellent  protection  against 
adenovirus  and  influenza  while  GAS  is  controlled  with  antibi¬ 
otic  treatment  and/or  prophylaxis.’  When  sulfonamide  antibiot¬ 
ics  were  discovered  the  US  military  employed  them  to  combat 
GAS  pharyngitis  and  the  severe  pneumonias,  rheumatic  fever, 
necrotizing  fasciitis  and  other  morbidity  that  can  result  from 
invasive  GAS  infection.^  This  practice  was  initially  effective  but 
within  1  year  resistant  strains  evolved  and  spread,  rendering  sul¬ 
fonamides  ineffective.^  When  penicillin  was  discovered,  it  too 
was  deployed  to  decrease  the  impact  of  GAS.  Penicillin  was  used 
prophylactically  on  entire  classes  of  trainees  in  order  to  control 


both  active  infections  and  asymptomatic  colonizations  that  act 
as  reservoirs. 

The  practice  of  active  prophylaxis  has  continued  for  over 
50  years  and  remains  effect ive.’  "**^  Efforts  to  limit  or  suspend 
prophylaxis  of  US  military  recruits  for  periods  of  more  than  a 
few  months  have  often  resulted  in  GAS-associated  outbreaks  of 
pharyngitis  and  ARD.^'^’*^  Penicillin  resistance  has  never  arisen 
in  GAS,  which  may  be  adaptively  incapable  of  acquiring  and/ 
or  expressing  the  most  common  penicillin  resistance  factors.’^'^ 
Between  7  and  30%  of  recruits  are  not  treated  with  penicillin 
due  to  demonstrated  or  potential  allergies.  Some  recruit  facilities 
give  these  individuals  other  antibiotics,  usually  macrolides  such 
as  erythromycin  or  azithromycin.’  ’^’  While  untreated  subp^opula- 
tions  experience  occasional  outbreaks  of  GAS,  they  are  partially 
protected  by  prophylaxis  of  the  majority  of  the  population.’  ’”” 
This  dynamic  is  similar  to  that  of  herd  immunity  (transmission 
interference  by  the  proportion  of  individuals  rendered  nonsus- 
ceptible  through  preexisting  immunity). 
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GAS  strains  are  highly  variable  in  their  rate  of  transmission, 
antibiotic  resistance,  range  of  disease  and  specific  antigenic¬ 
ity  3.6.7.12-15  is  most  commonly  characteri/xd  by  serotype  and 
by  antibiotic  resistance  profile.  The  most  common  method  of 
direct  serological  identification  is  M  typing,’^  but  this  has  been 
widely  replaced  by  inference  of  M  type  from  emm  gene  sequence 
analysis,’"'  or  more  indirectly  inferred  by  various  forms  of  multi¬ 
locus  sequence  typing.”’  All  of  these  methods  have  shown  broad 
concordance  to  the  degree  that  they  overlap  in  coverage. 

Analysis  of  the  sample  set  collected  for  this  study  showed  very 
strong  correlations  between  e?n?n  type  and  antibiotic  resistance 
phenotypes,  as  well  as  strong  tendencies  lor  specific  types  to  be 
aSxSOciated  with  outbreaks.”^  Previous  studies  report  similar  asso¬ 
ciations,  supporting  the  temporal  and  spatial  stability  of  observed 
correlations.  The  strong  correlation  between  serotype  and  phe¬ 
notype  suggests  possible  clinical  uses  of  rapid  genotyping  for  the 
inference  of  antibiotic  susceptibility  and  the  prediction  of  epide¬ 
miological  and  clinical  severity.”^ 

Turnover  between  different  serotypes  is  likely  driven  by  cycles 
of  serotype-specific  herd  immunity  in  the  general  population 
that  tend  to  favor  recently  rare  strains.  A  larger  proportion  of  the 
population  is  susceptible  to  recently  rare  types,  arfd  hence  they 
are  more  readily  transmitted.  Dynamic  epidemiology,  character¬ 
ized  by  rapid  turnover  between  dominant  strains,  is  the  recog¬ 
nized  pattern  for  GAS.'”"'’  However,  the  scale  on  which  these 
turnovers  occur  has  not  been  well  studied.  Complete  turnover 
between  single  dominant  strains  was  seen  in  a  small,  somewhat 
isolated  population  of  children,'’  while  less  succinct  dynamic 
events  were  also  seen  in  military  recruits.'”  The  authors  of  those 
studies  attributed  the  turnover  events  to  different  factors  herd 
immunity  and  close  contact  in  the  first  case,  constant  turnover  of 
the  population  in  the  second.  Here,  we  observe  multiple  sites  over 
a  period  of  6  years  to  address  the  potential  connection  of  turnover 
events  on  a  larger  scale  and  the  impact  of  these  events  on  local 
changes  in  antibiotic  resistance  and  virulence. 

We  present  6  years  of  GAS  surveillance  data  including  site- 
specific  strain  distribution  time  series,  GAS  disea.se  rate  data  and 
temporal  patterns  of  macrolide  resistance  from  ten  US  military 
recruit  training  facilities.  We  explore  the  relationship  between 
strain  distribution  patterns  at  multiple  sites.  Finally,  the  contri¬ 
butions  of  time,  location,  prophylaxis  regimen  and  emm  type  to 
macrolide  resistance  patterns  are  discu.ssed. 

Results 

Molecular  epidemiology  of  GAS  among  military  recruits. 
Temporal  .serotype  di.stributions  and  a.ssociated  rates  of  ARD  and 
macrolide  resistance  are  shown  in  Figures  1  and  2,  and  the  leg¬ 
end  defining  serotype  identities  and  overall  distribution  is  shown 
in  Figure  3.  The  as.sociated  patterns  of  antibiotic  use  at  the  stud¬ 
ied  sites  arc  shown  in  Figure  4. 

The  simultaneity  of  specific  .serotype  emergences  and  turn¬ 
overs  is  clear  from  visual  examination.  From  2002  through  2004, 
M75,  M6  and  M3  were  dominant  at  the  two  individual  sites  and 
in  the  composite  group.  These  strains  became  universally  less 
common  from  2005  through  2007.  They  were  replaced  by  M44 


and  M5  at  FLW  and  MCRD-Pl,  by  Ml  18  at  all  sites,  and  by  MI8 
at  FLW  and  in  the  composite  group. 

In  order  to  address  the  significance  of  these  ob.served  correla¬ 
tions  we  binned  data  across  months  and  years  and  tested  each 
emm  type  for  temporal  correlation  of  dominance  patterns  among 
the  two  primary  sites  and  the  composite  group  (see  Table  1). 
Positive  correlations  were  common  and  often  significant  for  well- 
sampled  emm  types,  while  negative  correlations  were  invariably 
nonsignificant,  supporting  our  subjective  visual  interpretation. 

Strain  turnover  and  antibiotic  resistance.  Only  a  few  of  the 
many  emm  types  of  GAS  are  responsible  for  the  vast  majority 
of  antibiotic  re.si.stancc.  M75  was  responsible  for  67%  of  the 
macrolide  resistance  observed  during  the  study  period  despite 
representing  only  10%  of  the  isolates  tested  for  resistance’'’ 
and  macrolide  resistance  did  not  show  any  significant  tempo¬ 
ral  and  spatial  trends  independent  of  the  distribution  of  M75 
(see  Figs.  1  and  2).  The  specific  associations  between  emm  type 
and  antibiotic  re.sistance  for  this  sample  set  has  been  previously 
published. 

Strain  turnover,  virulence  and  rate  of  GAS  pharyngitis. 
Rate  data  are  .superirnpo.sed  on  the  serotype  distribution  graphs 
(Fig.  1)  for  time  periods  during  which  both  site  population  and 
GAS-a.s.sociated  disease  rate  were  measured.  In  late  2005  and  2006, 
there  was  a  marked  increase  in  GAS  morbidity  among  recruits. 
The  noted  outbreaks  of  invasive  di.sease  and  ARD  in  2005-2006 
were  closely  associated  with  the  dominance  of  M5-  The  emer¬ 
gence  of  this  apparently  virulent  and  highly  transmissible  strain, 
which  was  associated  with  many  cases  of  invasive  di.sease  and  at 
least  two  fatalities,  preceded  a  temporary  shortage  of  Bicillin  that 
cau.sed  a  suspension  of  prophylaxis  at  many  sites  in  late  2006. 

Discussion 

Molecular  epidemiology  of  GAS  among  military  recruits. 
The  data  reveal  highly  dynamic  epidemiology,  albeit  against  a 
more  diverse  background  than  .seen  in  smaller  and  more  isolated 
populations.'’  While  it  has  been  suggested  that  turnovers  in 
military  populations  simply  represent  the  effect  of  training  class 
turnover  events,'”  the  data  suggest  otherwise.  The  sites  studied 
here  experience  constant  training  class  turnover,  with  approxi¬ 
mately  10%  of  the  recruit  population  being  replaced  weekly, 
yet  specific  emm  types  often  dominate  specific  sites  for  periods 
of  several  months.  Different  sites  often  retain  different  distri¬ 
butions  oi  emm  types  for  months  at  a  time,  despite  constantly 
receiving  recruits  from  the  .same  general  population,  suggesting 
that  transmijision  of  specific  strains  is  maintained  within  sites 
for  periods  of  several  months.  For  example,  M75  is  dominant 
at  MCRD-PI  while  M3  is  dominant  at  FLW  from  February 
through  May  2004  (.see  Fig.  1). 

On  the  other  hand,  these  two  geographically  distinct  training 
.sites  (approximately  500  miles  apart)  experience  many  dynamic 
epidemiological  events  in  parallel.  While  locally  endemic  strains 
may  be  transmitted  within  sites  on  a  monthly  scale,  turnovers 
on  a  longer  (yearly)  .scale  are  likely  reflective  of  widespread  .strain 
turnover  events  that  extend  across  multiple  military  training  sites 
and  po.ssibly  the  entire  general  population. 


248 


Virulence 


Volume  1  Issue  4 


Figure  1.  emm  type  distribution,  GAS  disease  rate  and  erythromycin  resistance  rate  over  6  years  at  two  sites.  Temporally  aligned  emm  type  distribu¬ 
tions  for  two  sites  that  consistently  contributed  samples.  Full  names  and  locations  of  all  sites  are  given  in  Acknowledgments,  emm  type  color  key 
appears  in  Figure  3. 


Ft  Benning,  Ft  Sill.  Ft  Gordon,  Great  Lakes.  Ft  Knox,  LL  AFB,  MCRD  SD  and  Madigan 

(Composite  Distribution) 


Figure  2.  Composite  emm  type  distribution  and  erythromycin  resistance  rate  over  6  years  at  eight  sites.  Temporally  aligned  emm  type  distribution  for 
a  composite  of  eight  sites  which  contributed  samples  sporadically.  Full  names  and  locations  of  all  sites  are  given  in  Acknowledgments,  emm  type  color 
key  appears  in  Figure  3. 
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Figure  3.  Overall  emm  type  distribution  at  ten  sites  and  key  to  Figures  1  and  2. 
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Figure  4.  Prophylaxis  regimens  used  to  prevent  and  control  GAS  outbreaks.  Automatic  accession  dose  =  prophylaxis  upon  entry  to  training.  Multidose 
schedules  may  vary  by  site,  generally  every  3  or  4  weeks.  Surveillance  indicates  sporadic  prophylaxis  based  on  observes  GAS  disease  rates.  Site  names 
are  color  coded  by  service  (see  Acknowledgments). 
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Strain  turnover  and  antibiotic  resistance.  The  data  support 
the  hypothesis  that  resistance  dynamics  result  from  natural  fluc¬ 
tuations  in  the  relative  dominance  of  serotypes  that  differ  in  terms 
of  their  inherent  (and  relatively  invariable)  resistance  phenotypes, 
pigores  1  and  2  show  this  clearly,  as  macrolide  resistance  can 

seen  to  cycle  with  the  dominance  of  M75  at  all  studied  sites. 
This  is  the  opposite  of  what  is  seen  in  many  bacteria,  such  as 
Streptococcus  pneumoniae, and  provides  evidence  for  relatively 
strong  genetic  stability  (linkage)  in  GAS. 

FLW  used  year-round  penicillin  prophylaxis  for  all  nonallergic 
recruits  but  never  used  macrolides  (or  other  antibiotics)  for  pro¬ 
phylaxis  of  penicillin-allergic  recruits  during  the  period  described 
in  this  study.  From  1998  through  early  2006,  MCRD-Pl  used 
prophylactic  penicillin  during  the  peak  GAS  season  (October 
1  through  March  31),  and  during  those  periods  they  also  used 
oral  macrolides  (azithromycin)  for  penicillin-allergic  recruits. 
Both  sites  experienced  dominance  of  the  primary  macrolide- 
resistant  strain  M75  in  2002,  independent  of  the  fact  that  one 
site  was  actively  using  macrolides  while  the  other  was  not.  M75 
did  remain  dominant  longer  at  MCRD-PI  than  at  FLW,  but  use 
of  macrolide  prophylaxis  did  not  prevent  the  subsequent  replace¬ 
ment  of  M75  with  macrolide-susceptible  strains  such  as  M5, 
Ml  18  and  M44  (see  Fig.  1). 

From  October  1,  2006,  through  March  31,  2007,  MCRD-Pl 
used  azithromycin  prophylaxis  for  all  recruits,  due  to  a  national 
shortage  of  the  Bicillin  (penicillin  G  benzathine)  formulation 
normally  used  for  this  purpose.  This  active  use  of  macrolides  was 
not  followed  by  the  emergence  or  dominance  of  macrolide-resis- 
tant  strains,  and  was  not  associated  with  the  appearance  of  new 
resistance  in  previously  sensitive  serotypes. 

The  pattern  of  specific  antibiotic  resistance  phenotypes  being 
associated  with  specific  emm  types  held  generally  true  for  all  emm 
types  and  all  antibiotics, with  one  exception.  The  only  hospital 
surveyed  in  this  study  (Madigan  Army  Medical  Center)  contrib¬ 
uted  just  five  isolates,  yet  these  revealed  very  high  rates  of  eryth¬ 
romycin,  tetracycline  and  clindamycin  resistance  (see  Table  2). 
Furthermore,  these  isolates  were  of  four  different  emm  types,  and 
these  emm  types  did  not  generally  display  the  same  resistance 
patterns  in  the  rest  of  the  sample  set.  This  suggests  genetically 
mobile  resistance  factors  and  active  selection  for  resistance  are 
playing  a  role  in  the  hospital.  These  phenomena  were  not  evi¬ 
denced  in  nonhospital  environments’^^  (see  also  Figs.  1  and  2), 

Strain  turnover,  virulence  and  rate  of  GAS  pharyngitis.  The 
data  are  consistent  with  previous  reports  suggesting  that  only  a 
small  number  of  the  known  emm  types  are  generally  associated 
with  high  virulence  and  invasive  characteristics,  with  Ml,  M3 
and  M5  foremost  among 

In  much  the  same  way  that  the  rise  and  fall  of  M75  appears 
responsible  for  recent  spikes  in  macrolide  resistance,  the  data 
suggest  that  the  observed  spikes  in  GAS-associated  ARD  and 
invasive  disease  were  associated  with  the  shifting  dominance  of 
a  particularly  virulent  strain,  M5.  The  association  between  this 
serotype  and  the  widespread  rate  increases  seen  among  trainees 
from  2002  through  2007  supports  the  general  hypothesis  that 
virulence  is  strongly  linked  to  serotype.  The  fact  that  M5,  which 
did  not  display  any  antibiotic  resistance  phenotypes,  was  able  to 


Table  1.  Significance  of  temporal  correlation  of  emm  type  frequency 


between  sites 

MCRD  Parris  Island 

Fort  Leonard  Wood 

By  month 

By  month 

Fort  Leonard  Wood 

5  (38%;  0.002) 

b 

75  (31%;  0.01) 

b 

44  (58%;  <0.0001) 

b 

Composite 

12  (40%;  0.001) 

18(62%;  <0.0001) 

By  year 

1  (34%;  0.004) 

By  year 

Fort  Leonard  Wood 

44  (97%;  0.002) 

b 

Composite 

3  (97%;  0.02) 

18  (99%;  0.0003) 

118(96%;  0.002) 

6(91%;  0.01) 

Data  shown  as  "emm  type  (correlation  %;  Pearson  p  value)".  All 
significant  correlations  were  positive  (distributions  changed  more 
synchronously  than  expected  by  chance).  Nonsignificant  correlation 
values  and  associated  p  values  are  not  shown.  All  emm  types  were  ana¬ 
lyzed  that  represented  at  least  3%  of  the  total  Isolates  (those  shown  in 
color  in  Fig.  3).  The  most  heavily  sampled  emm  types  offered  the  most 
positive  correlations,  suggesting  that  significance  may  be  limited  by 
sample  size  for  other  emm  types.  ‘’Blank  cells  are  self  comparative. 

cause  outbreaks  at  a  site  (MCRD-PI)  where  70%  of  the  recruits 
were  receiving  Bicillin  prophylaxis  twice  during  training  demon¬ 
strates  the  ongoing  threat  of  GAS  to  US  military  trainees.  The 
diverse  antibiotic  resistance  and  virulence  profiles  of  different 
emm  types  necessitate  flexible  and  responsive  prophylaxis  pro¬ 
grams,  and  demonstrate  the  potential  value  of  rapid  genotyping 
as  a  tool  for  the  prediction  of  clinical  and  epidemiological  severity 
and  for  informing  treatment  and  prophylaxis  decisions. 

Materials  and  Methods 

Diagnostic  samples  were  collected  by  culture  on  blood  agar  plates 
from  patients  with  pharyngitis  or  other  apparent  GAS  disease. 
Collection  was  performed  by  on-site  medical  staff  as  part  of  rou¬ 
tine  diagnostic  procedures  at  10  US  military  facilities  between 
January  2002  and  December  2007  (collection  personnel  were  not 
members  of  the  author  group).  Cultured  isolates  remaining  from 
diagnostic  procedures  were  sent  to  the  Naval  Health  Research 
Center  with  minimal  demographic  information  for  purposes 
of  strain  typing  and  antibiotic  resistance  testing  under  Naval 
Health  Research  Center  Institutional  Review  Board  protocol 
NHRC. 2001. 0008.  Consent  was  waived  because  the  work  was 
done  entirely  on  deidentified,  already  existing  cultured  isolates 
collected  in  the  process  of  routing  diagnostic  testing,  and  results 
were  not  used  for  diagnosis  or  patient  management.  This  research 
has  been  conducted  in  compliance  with  all  applicable  federal  reg¬ 
ulations  governing  the  protection  of  human  subjects  in  research. 

Serotyping  data  from  the  preceding  5  years  of  this  program, 
covering  many  of  the  same  sites,  have  been  described  else¬ 
where."  Collection  was  irregular,  being  dependent  on  available 
time,  awareness  of  and  access  to  apparent  cases  of  GAS  disease, 
and  the  presence  or  absence  of  such  disease  at  any  given  time. 
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Table  2.  Temporal  and  geographic  distribution  of  antibiotic  resistance  phenotypes  among  group  A  streptococcus  isolates 


Erythromycin 

Chloramphenicol 

Ofloxacin 

Tetracycline 

Year 

N 

% 

N 

% 

N 

% 

N 

% 

2002 

103 

16% 

80 

3% 

80 

15% 

103 

3% 

2003 

171 

24% 

171 

9% 

171 

15% 

171 

7% 

2004 

77 

25% 

77 

5% 

76 

3% 

77 

7% 

2005 

70 

4% 

70 

19% 

70 

3% 

70 

4% 

2006 

271 

3% 

63 

14% 

63 

0% 

271 

1% 

2007 

110 

7% 

0 

0 

110 

2% 

Site 

N 

% 

N 

% 

N 

% 

N 

% 

FLW 

274 

8% 

152 

11% 

151 

7% 

274 

3% 

Ft.  Benning 

31 

13% 

31 

3% 

31 

13% 

31 

23% 

Ft.  Sill 

85 

7% 

84 

13% 

84 

11% 

85 

4% 

Ft.  Gordon 

21 

0% 

0 

0 

21 

0% 

Great  Lakes 

92 

7% 

37 

11% 

37 

11% 

92 

2% 

Ft.  Knox 

39 

0% 

0 

0 

39 

0% 

Lackland 

44 

0% 

37 

3% 

37 

19% 

44 

5% 

MCRD-Pl 

140 

34% 

100 

y  10% 

100 

6% 

140 

1% 

MCRD-SD 

71 

7% 

20 

/  0% 

20 

5% 

71 

6% 

Madigan 

5 

80% 

0 

0 

5 

20% 

Clindamycin 


N 

103 

171 

77 

70 
271 

no 

N 

274 

31 

85 

21 

92 

39 

44 

140 

71 
5 


% 

1% 

2% 

0% 

0% 

1% 

6% 

% 

0% 

0% 

0% 

0% 

1% 

0% 

0% 

3% 

3% 

100% 


Distribution  is  shown 
includes  both 
section. 


shown  for  all  antibiotics  to  wnicn  uicr  luicii  - -  ,  .  .  .  *  ,  •  ^ _ _ 

intermediate  and  resistant  phenotypes  per  CLSI  definitions.  Full  names  and  locations  of  all  sites  are  given  in  the  Acknowledgme 


Two  facilities,  Fort  Leonard  Wood,  Missouri  (FLW)  and  Marine 
Corps  Recruit  Depot  Parris  Island,  South  Carolina  (MCRD-Pl) 
consistently  provided  samples  over  most  of  the  study  period, 
while  other  sites  submitted  samples  sporadically.  FLW  consis¬ 
tently  collected  both  recruit  population  size  (denominator)  data 
and  GAS-associated  illness  (numerator)  data,  as  did  MCRD-Pl 
from  late  2003  through  early  2004  and  from  2005  through  2007. 
This  data  allowed  calculation  of  GAS-associated  illness  rates,  as 
depicted  in  Figures  1  and  2. 

Positive  cultures  were  resuspended  in  tryptic  soy  broth  with 
13%  glycerol,  shipped  frozen  on  dry  ice  to  the  Naval  Health 
Research  Center  (NHRC),  and  stored  at  -70"C.  Traditional  emm 
sequencing  was  employed  to  type  387  samples.  The  emm  types 
of  the  remaining  415  isolates  were  inferred  by  PCR  electrospray 
ionization  mass  spectrometry  (PCR/ESI-MS)  analysis  (a  rapid 
form  of  multilocus  genotyping'**)  with  reference  to  correlations 
between  PCR/ESl-MS  identity  and  emm  sequence  identity  at 
the  same  site  during  the  same  year.’^^  Antibiotic  resistance  profiles 
were  generated  using  standard  culture-based  inhibition  assays  as 
previously  described. 

In  an  effort  to  discern  the  geographic  scale  of  observed  strain 
turnover  events,  the  temporal  patterns  of  serotype  distribution  at 
the  two  consistently  sampled  sites  were  compared  with  each  other 
and  with  a  third  composite  group  representing  less  densely  sam¬ 
pled  sites.  Pearson  correlation  analyses  were  applied  to  measure 
the  statistical  significance  of  temporal  correlations  in  dominance 
of  individual  emm  types  between  pairs  of  sites  on  monthly  and 
yearly  scales. 


General  Conclusions 

Strain  distributions  appear  to  shift  in  a  correlated  fashion 
across  multiple  sites.  This  suggests  that  the  observed  dynamic 
epidemiology  is  reflective  of  strain  turnover  events  occurring 
on  larger  population  scales.  Macrolide  resistance  during  the 
study  period  came  primarily  from  M73,  and  rates  of  mac¬ 
rolide  resistance  were  closely  correlated  with  the  temporary 
dominance  patterns  of  M75.  The  same  appears  to  be  true 
of  tetracycline  resistance,  with  almost  all  resistance  com¬ 
ing  from  3  emm  types:  M38,  M68  and  M77’‘'  (temporal  data 
for  tetracycline  resistance  not  shown).  Similarly,  high  rates 
of  GAS-associated  ARD  and  invasive  disease  were  strongly 
associated  with  M3  and  M3,  both  of  which  have  previously 
been  associated  with  high  virulence  and  outbreak  potential. 
We  therefore  argue  that  changes  in  GAS  morbidity  and  mor¬ 
tality,  as  well  as  changes  in  antibiotic  resistance  rates,  are  pri- 
marily  the  re.sult  of  geographically  widespread  changes  in  the 
relative  dominance  of  a  few  types  with  unique  phenotypes. 
Antibiotic  resistance  does  not  seem  to  he  strongly  affected 
by  prophylaxis  regimens.  Based  on  this,  we  argue  that  emm 
type-specific  surveillance  could  offer  a  high  degree  of  predic¬ 
tive  value  related  to  clinical  impact,  epidemiological  sever¬ 
ity  and  antibiotic  resistance.  Tracking  resistance  patterns  for 
the  purposes  of  controlling  and  understanding  the  spread  o 
resistance  should  also  be  done  on  an  emm  type-specific  basis 
or  any  significant  trends  will  be  obscured  by  the  indirect 
effects  of  strain  turnover  dynamics. 
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